Local phase space control and interplay of classical and quantum effects in 
dissociation of a driven Morse oscillator 
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This work explores the possibility of controlling the dissociation of a monochromatically driven 
one-dimensional Morse oscillator by recreating barriers, in the form of invariant tori with irrational 
winding ratios, at specific locations in the phase space. The control algorithm proposed by Huang et 
al. (Phys. Rev. A 74, 053408 (2006)) is used to obtain an analytic expression for the control field. 
We show that the control term, approximated as an additional weaker field, is efficient in recreating 
the desired tori and suppresses the classical as well as the quantum dissociation. However, in the case 
when the field frequency is tuned close to a two-photon resonance the local barriers are not effective 
in suppressing the dissociation. We establish that in the on-resonant case quantum dissociation 
primarily occurs via resonance-assisted tunneling and controlling the quantum dynamics requires a 
local perturbation of the specific nonlinear resonance in the underlying phase space. 



I. INTRODUCTION 

For over three decades the one-dimensional driven 
Morse oscillator [lj has served as a fundamental model 
to understand and elucidate the dissociation mechanism 
of diatomic molecules. The continued interest in this, 
seemingly simple, system is due to two main reasons. 
First, the hope is that insights into the mechanism can 
be utilized to understand infrared multiphoton dissocia- 
tion of polyatomic molecules [2, U 13 an d related phenom- 
ena including vibrational predissociationQ and mode- 
specific dynamics [1, 0|- Second, at present the focus of 
researchers is increasingly shifting from gaining mech- 
anistic insights to controlling 0, 0, [l(J, [ijj the various 
processes and in this regard a firm understanding of the 
underlying mechanisms is essential. Therefore, it is not 
entirely surprising that the driven Morse system has been 
studied in great detail from the quantum, classical, and 
semiclassical[l2j perspectives and with an equally diverse 
choice for the field - monochromatic (l3l [Til flU Ha . Il7j . 18, 



211 1221 . bichromatic (with relative phase) [23t 124 
^j3j2<| [U, chirpedHH, [H, Si El, shaped 
pulses [2 ll I261 . l3o| . and stochastic noise[3a]. More re- 
cently, the dynamics of a Morse oscillator under the influ- 
ence of external fields has become relevant in the context 
of models for quantum computing based on molecular 
vibrations [3 7j. 

A majority of the studies have addressed the prob- 
lem from a classical-quantum correspondence viewpoint; 
a knowledge of the regimes where classical or quantum 
mechanisms are appropriate and regimes where they co- 
exist and compete is crucial for control[38j. Several im- 
portant insights have originated from such efforts which 
have established that molecular dissociation, in analogy 
to multiphoton ionization of atoms [39l l4fj| . [4l| . occurs 
due to the system gaining energy by diffusing through 
the chaotic regions of the phase space. For example, an 
important experimental study by Dietrich and Corkum 
has shown [42j|. amongst other things, the validity of the 
chaotic dissociation mechanism. Thus, the formation of 
the chaotic regions due to the overlap [43j of nonlinear 



resonances (field- matter) , hierarchical structures [441 l45j 
near the regular-chaotic borders acting as partial barri- 
ers, and their effects on quantum trans port 1461 have been 
studied in a series of elegant papers fiH l47l. |48| . A general 
consensus, atleast for the one-dimensional driven Morse 
system, is that classical-quantum correspondence holds 
up rather well except in the regimes of quantum multi- 
photon resonances [H, [H, [I?], O H3| . 

A recurring theme in many of the works on driven 
Morse system has to do with enhancing the dissoci- 
ation. The search for ways to efficiently dissociate 
the molecule has led to a variety of suggestions like 
bichromatic fields with the relative phase as a control 
knob[H [H [H, [H [3(1 frequency-chirped fields [ij , and 
resonant stimulation[49|| . However, there are instances 
wherein one is interested in suppressing the dissociation 
rather than enhancing it. This is important, for exam- 
ple, in the context of vibrational quantum computing (37| 
where loss of population into states other than the states 
of interest compromises the efficiency of the quantum 
gates. Another example comes from coupled Morse os- 
cillator systems where it might be neccessary to keep 
one of the modes 'quiet' in order to carry out mode- 
specific dynamics A powerful approach to implement 
such constraints on the system comes from optimal con- 
trol theory (Hoj] (OCT) and indeed driven Morse oscil- 
lator systems provide an ideal testbed for OCT-based 
schemes (51). Yet, in our opinion, it is worthwhile ad- 
dressing the issue from a classical-quantum correspon- 
dence perspective as well. Not only is it natural, given 
the extensive insights that classical mechanics can pro- 
vide, but it might also provide a useful way of decoding 
information buried in an otherwise complicated optimal 
field coming out of an OCT calculation. Similar con- 
siderations are at the heart of several works [52| aimed at 
understanding the dynamical origins of the control fields. 

Since a detailed understanding of the role of various 
phase space structures in the driven Morse system al- 
ready exists, is it possible to use the phase space infor- 
mation to control the dissociation using additional, hope- 
fully simple, fields? Recently[53], a similar question was 
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addressed by Huang et al. in the context of suppress- 
ing the multiphoton ionization of atomic systems. Us- 
ing methods [54I developed in a different context, it was 
found that the ionization process could be suppressed by 
rebuilding some of the broken invariant tori at carefully 
chosen locations in the phase space. Inspired by their 
approach, and noting the mechanistic similarities be- 
tween molecular dissociation and atomic ionization(55j. 
in this work we attempt to control the dissociation of 
a monochromatically driven Morse oscillator using the 
local control algorithm. In their study, Huang et al. 
focused [53fl only on the classical aspects of suppressing 
the ionization. It is, however, important to ask if the 
classical barriers are effective quantum mechanically as 
well since it is not immediately clear that local barriers in 
the phase space translate to local suppression of quantum 
dynamics. We address this issue using the driven Morse 
system and show that phase space barriers, especially 
cantori, do inhibit both classical and quantum dissocia- 
tion. As one would expect, such good classical-quantum 
correspondence fails in the case of two-photon resonance. 
However, we show that the complication comes from a 
subtle interplay between classical and quantum mecha- 
nisms with resonance-assisted tunneling [5^. l57l. [58j play- 
ing a key role. 

We begin by describing some of the salient features of 
the driven Morse oscillator in section [II] After a brief 
description of the methodology, section IIIII contrasts the 
dissociation dynamics in the off-resonant and on-resonant 
situations and a specific initial Morse state is identified to 
be subjected to the local control strategy. In section HVl 
we give a brief summary of the local control method re- 
sulting in an analytic form of the control field. A simpli- 
fied control field, appropriate for classical-quantum cor- 
respondence studies, is obtained. The efficiency of the 
simplified control term in recreating various cantori barri- 
ers in phase space and hence controlling the classical and 
quantum dissocation dynamics is shown and discussed in 
section [Vj In the same section we illustrate the impor- 
tance of resonance- assisted tunneling in the on-resonance 
regime. Finally we conclude in section I VII with some 
comments on the method, possible generalizations, and 
future applications. 



II. MODEL HAMILTONIAN 

The driven Morse oscillator, modelling the dissociation 
of a diatomic molecule by linearly polarized laser fields, 
is described[l5| by the Hamiltonian 

H(x,p;t) = H (x,p) - \ 1 n(x)cos(u> F t), (1) 

with the unperturbed Hamiltonian 



Ha(x,p) = ^p 2 + D [l 



-Qt(x— %e)] 



(2) 



corresponding to a one-dimensional Morse oscillator. It 
is well known that H Q provides a good model for describ- 
ing the anharmonic vibrations of diatomic molecules with 



D 0l a, x e , and M being the dissociation energy, range of 
the potential, equilibrium position, and the reduced mass 
of the molecule respectively. The bound eigenstates and 
eigenvalues corresponding to the Hamiltonian Hq(x,p) 
can be expressed, with z = 2ae~ a ( x ~ x '\ as[59| 
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(3b) 



where L 2bv (z) is the generalized Laguerre polynomial, 
a = y/2MD a /ah and, b v = ^-2E v M/ah. 

A fit [13] to the ab-initio data on HF yields the following 
form for the dipolc 



(i(x) = Axe 0X , 



(4) 



with A = 0.4541 a.u. and (5 = 0.0064 a.u. However, the 
linear approximation 



fi(x) 



(X - X e ) 
/J,{x e ) +dl(x- X e ), 



(5) 



with d\ « 0.33 a.u. is used in the current work since, 
for the moderate field intensity of interest, the qualita- 
tive nature of classical and quantum dynamics are unal- 
tered as compared to working with the dipole function 
in Eq. 21 Moreover, as seen later, the linear form allows 
for a relatively easier implementation of the local control 
algorithm [54| in terms of deriving analytic expressions 
for the control field. 

Given that this work focuses on suppressing disso- 
ciation by creating robust Kolmogorov-ArnoPd-Moser 
(KAM) tori in the phase space, the action-angle variables 
(J, 6) of the unperturbed Morse oscillator 



J 



2MD 



1 



-sgn(p) cos 1 



l-E 



a(x—x e ) 



E 



(6a) 
(6b) 



are a convenient and natural representation to work with. 
In terms of (J, 9) the Hamiltonian in Eq.[T]can be written 
down as 



H(J,9;t) = H (J) -ev(J,9;t), 
with e = Ai<ii and 



H (J) = w [J 
v(J,0;t) = 2 



4D 

OO 

Va(J) + ^V n (J) cos(n0) 



(7) 



(8) 



n=l 



cos(wp£), 



being the zeroth-order Hamiltonian and the perturbation 
respectively. In the above equations, E — Hq/Dq < 1 



denotes the dimensionless bound state energy, loq — 
[2a 2 Dq/M) 1 ' 2 is the harmonic frequency, and sgn(p) = 
1 for p > 0, sgn{p) = —1 for p < 0. The Fourier coef- 
ficients Vq(J) and V n (J) are known analytically [HI [47| 
and given by 



Vo(J) 



V n (J) 



4a 
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D + V A) - D E(J) 

2(D -E(J)) 



(-1) 



n+1 



on 



V D oE(J) 



Do + V D o - D oE(J) 



(9a) 



(9b) 



Note that the classical nonlinear frequency of the Morse 
oscillator is given by 



U0{J) = -7TT = WO ( 1 - 7T7^~J 



dJ 



2D 



(10) 



Throughout this work, we use atomic units for 
the various parameters including the field parameters 
(di, Ai, Lop) and measure time in units of the field period 
t = 2tt/luf. In particular, we choose D — 0.225, a = 
1.174, x e = 1.7329 and, M = 1744.59 corresponding to 
the Hydrogen Fluoride (HF) molecule [ToT]. The poten- 
tial well supports a total of Nb = 24 bound states. The 
laser field amplitude is fixed at Ai = 0.0287 a.u. (~ 30 
TW/cm 2 ), and thus we are studying dissociation under 
a moderate intensity field. 
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FIG. 1: (Color online) Classical dissociation probabilities for 
the Morse oscillator states v — 10 (circles), 11 (squares), 12 
(triangles) , and 13 (diamonds) with driving field frequency (a) 
lof = 0.0129 (b) 0.0178 and fixed field amplitude A = 0.0287 
(~ 30 TW/cm 2 ). The corresponding quantum results are 
shown in the right column as (c) and (d) respectively. In case 
of uif — 0.0178 an additional state v = 14 (line, no symbol) is 
also shown. Note that in this figure and all subsequent figures 
the various parameter values are in atomic units. 



III. CLASSICAL AND QUANTUM 
DISSOCIATION DYNAMICS 

Although one can choose different classes of initial 
states for the study, in this work the initial states are 
chosen to be the zeroth-order Morse eigenstates x„ given 
by Eq. I3al The initial states are time evolved on a grid 
using the well established split-operator method[61| in- 
volving the short-time propagator 



U(At) = exp 



2h 



exp 



At f, 



exp 



.My 

2ft 



(11) 

with T and V denoting the kinetic and potential en- 
ergy operators respectively. The time-step was set to 
At = 5 x 10 _3 r to ensure convergence of the dissociation 
probabilities over the timescales of interest of about 500r. 
As is usual, unphysical reflection at the grid boundaries 
is avoided by employing an optical potential pfil. |62| 



V opt (x) 



iV 



(1 + el-^-**)/'']) ' 



(12) 



with parameters V = 0.02, r) = 0.35, and x* = 16.74. 
The introduction of V op t smoothly damps the outgoing 
wavefunction and does not modify the time-evolution of 
the bound states. The quantum dissociation probability 



is then calculated as 



(13) 



where Nb is the number of bound states. 

In order to compare and contrast the quantum disso- 
ciation dynamics with the classical dissociation dynam- 
ics we compute [63] the classical dissociation probabili- 
ties Pc(t) by choosing an ensemble of initial trajectories 
N to t with energy E v corresponding to the specific ini- 
tial Morse state with the angle uniformly distributed in 
[— 7r, 7r]. During the time evolution, a trajectory is con- 
sidered to be dissociated when the compensated energy 



E r = 



1 

2M 



P 



LOF 



- sin(u> ft) 



A>[1 



)i 



(14) 

exceeds the Morse dissociation energy Do- The number 
of dissociated trajectories N^is at a given time arc deter- 
mined from the above criteria and the resulting classical 
dissociation probability is the fraction 



P c d{t) 



(15) 



and P$) as a function of time 



In Fig. Q] we show P£> 
for some of the high-lying Morse eigenstates for two spe- 
cific driving field frequencies of u>f = 0.0129 and 0.0178 
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For the specific Morse parameters and field strength, the 
lower ujp represents an off-resonant situation whereas the 
higher field frequency corresponds to a two-photon reso- 
nant case. These cases, which will be used to highlight 
the results, are selected since they represent two limits 
in which classical-quantum correspondence either holds 
(off- resonant) or does not hold (on-resonant). A compar- 
ison of the classical (Fig. [IJa)) and quantum (Fig. [TJc)) 
results in the off-resonant case reveals that the dissoci- 
ation probabilities monotonically increase with increas- 
ing vibrational excitation. However, P|> is considerably 
smaller as compared to Pp. The reasons for this are 
well known and can be explained based on the clas- 
sical phase space shown in Fig. E^b). Extensive clas- 
sical stickiness [63l| around the initial action J = 10.5 
(corresponding to the quantum initial state v = 10) 
leads to the reduced Pfj for this state. At the same 
time the density variation in the chaotic regions of the 
phase space is symptomatic of the existence of partial 
barriers - in this case corresponding to a cantorus with 
uj F /n a {J) = I+7- 1 with 7= (V5+l)/2w 1.618 being 
the golden mean. Based on earlier works [Isl 147} . it is 
reasonable to surmise that the quantum dissociation is 
blocked by the cantorus. On the other hand, results for 
the on-resonant case shown in Figs. [l[b) and (d) indi- 
cate a fairly nontrivial behavior. The quantum dissoci- 
ation probabilities are non-monotonic with initial states 
v = 10,14, having nearly identical Pp, dissociating far 
more than the state v — 13. Quantum mechanically, res- 
onant two photon transition of state v = 10 to v = 14 
leads to direct coupling with the continuum and hence 
enhances the dissociation of state v = 10. The state- 
to-state transition probabilities indicate fl5j] . not shown 
here, Rabi cycling between the states v = 10, 12, and 
14. Insights into such behavior can also be gained by 
studying the classical phase space structures as seen in 
Fig. [2jc) . In this on-resonant case J = 10.5 is essen- 
tially located around the 1 + 7 -1 cantorus and a promi- 
nent u>f/CIq(J) — 2/1 nonlinear resonance is observed. 
Clearly, the 2:1 resonance is the classical analog of the 
quantum two-photon resonance and must be playing a 
crucial role in the observed dissociation dynamics [17]. In 
the subsequent sections we will highlight the classical- 
quantum correspondence for both the off-resonant and 
on-resonant cases. 

In order to illustrate the key features of this work we 
focus on the dissociation dynamics of the Morse state 
v = 10 for the above mentioned field frequencies. The 
analysis, however, can be performed for any initial state 
and our specific choice is inspired by the earlier work of 
Brown and Wvatt[l5|. Moreover, for the HF molecule, 
under moderate intensity fields, state v = 10 is a good 
choice to illustrate the interplay between classical and 
quantum dissociation mechanisms. Figure UK a) provides 
the comparison of quantum and classical dissociation 
probabilities for v = 10 interacting with a field with fixed 
intensity and for varying choice of the field frequencies 
ujp- The quantum distribution exhibits peaks at certain 
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FIG. 2: (Color online) (a) Quantum (squares) and classi- 
cal (circles) dissociation probability for the Morse eigenstate 
v = 10 of HF as a function of the field frequency ujf at the 
final time 500r with fixed Aj = 0.0287. Panels (b) and (c) 
show the stroboscopic surface of section for two representative 
cases considered in this work with luf — 0.0129 and 0.0178 
respectively. The dashed black lines in (b) and (c) are at the 
classical action value J = 10.5. The inset in (a) shows the 
effect of building the 1 + 7" 1 cantorus barrier on the clas- 
sical (filled circles) and quantum (filled squares) dissociation 
probabilities. 

frequencies corresponding to resonant multiphoton tran- 
sitions. The classical dissociation profile rises with u>f, 
broadens and dies out smoothly at higher frequencies due 
to the transition from trapping of trajectories in KAM 
tori at low frequencies to trapping inside the resonance 
island regions at higher frequencies. These observations 
are rather general and a detailed interpretation has been 
given earlier by Nicolaides and coworkers {63j. 

We now pose several questions in the context of Fig. [5] 
Is it possible to correlate the changes in classical phase 
space structures with varying uip with the quantum 
dissociation probabilities in both off-resonant and on- 
resonant cases? What is the role, if any, of the classical 
non-linear resonances in regulating the decay of quantum 
states? Finally, and the main focus of this work, can one 
control the classical and quantum dissociation dynamics 
by creating suitable local barriers in the classical phase 
space? For the present system the answer is in the af- 
firmative and, as a preview to the rest of the paper, the 
inset to Fig. [SJa) shows the suppression of classical and 
quantum dissociation by locally creating a cantorus with 
winding number u>f/Q.o(J) = 1 + 7 _1 - We now turn 
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to the issue of local phase space control of dissociation 
which, as seen later, provides answers to the first two 
questions posed above as well. 



IV. CONTROL BY REBUILDING A KAM 
TORUS 

The phase spaces shown in Figs.[3Jb) and (c), and the 
discussion in the previous section, suggest that if one can 
rebuild some of the irrational tori, such as the 1 + 7 _1 
cantorus, locally in the phase space then it ought to be 
possible to suppress the dissociation. Given the close 
parallels between the atomic ionization and the system 
of interest to us i.e., molecular dissociation, we employ 
the same technique j53j. based on classical perturbation 
theory, to obtain an analytic expression for the control 
field in case of the driven Morse oscillator. Since the 
technique has been described in considerable detail in the 
earlier works [53HH3|. in what follows we provide the main 
results which are of relevance in the present context. In 
addition, note that we use the notation of Huang et al. 
for convenience as well as uniformity. 



A. Methodology 

To start with, the nonautonomous Hamiltonian 
is mapped into an autonomous one by considering 
(t(mod27r) , E) as an additional angle-action pair. De- 
noting the action and angle variables by A = (J, E) and 
6 = (9, t) we can write the original driven system Hamil- 
tonian as 



H(A,0)=H (A)-eV(A,0). 



(16) 



Note that the unperturbed (e = 0) invariant tori la- 
beled by the action A correspond to the frequency oj = 
dHo/dA = (Q,,L)p). The aim is to rebuild a nonreso- 
nant torus A = ( Jo, 0), k • u> ^ with integer k, which 
has been destroyed due to the interaction with the field. 
This can be done by adding a small control term f(9) to 
H(A, 6) yielding the control Hamiltonian 

H c (A,0) = H(A,d) + f(0), 

f{6) = -H(A -d e Tb(0),d), (17) 

with b{8) = H(A Q , 9) = J2 k &ke 4k e and T being a linear 
operator defined by 



rb(e)= - 



— e lke . 
zk • U> 



(18) 



Referring to the phase spaces shown in Figs.[2]Jb) and 
(c) it is clear that the classical action corresponding to 
the quantum initial state v = 10 is located between the 
primary resonances u>p:£lo = 1:1 and 2:1. Thus, in our 
case the aim is to try and rebuild tori with irrational fre- 
quency ratios between the two resonances. In particular, 



the golden mean tori (1 + 7~", integer n) are of spe- 
cific interest in the driven Morse system[l5|, H3|- In the 
absence of external fields, such an invariant torus with 
frequency f2 r is located at J r = {uJq — fi r )(2U /^o 2 )- We 
shift the action J = J — J r to focus on the specific region 
of the phase space and expand the autonomous Hamilto- 
nian to second order (exact for the Morse oscillator) in 
J. Following the methodology outlined above the control 
term is obtained as 



/(M = ^r(.d e Tb(6)) 2 + e(2C + C x ), 



4/J 

where we have denoted 



(19) 



Cn 



E 

k=l 



^-V Qk (Jr)(dgTb) k (co S u; F t), 



* - E E 



{ -^v nk (j r )(d e rb) k \ 

71=1 U"=l ' J 

x [cos(n9 + uipt) + cos(n9 — wpt)] 



(20) 



with 



V nk (J r ) EE (^L Vn (J+Jr) 

\dJ k / ,7 =0 

OO 

Td e b(9) = d e Tb(9)=-eJ2nV n (J r ) (21) 



cos(n9 + LOpt) cos(n9 — uipt) 
(nfl r + up) (nf2 r — ujp) 



For moderate field intensities i.e., small e, one can work 
with the 0(e 2 ) approximation (leading order) to the con- 
trol term in Eq. [19] given by 



f a (9,t) = ^{doTbf -2eV 0X {d B Tb)cos{u>pt) 
- e((J,9;t)d e Tb(9), 



(22) 



and it can be shown that 

2 /2w, 



Voi - 







3 n {tuo-n^- 1 



v ( 1 \ n + 1 ( -^L—\ (° 



(23) 



OO 

C(J,8;t) = ^2 V nl (J r )[[cos(nd + UFt) + cos(n6 - w F t)]. 



It is important to note that the perturbative treatment 
carried out to derive the control field f(6, t) breaks down 
when nQ r i=s up. Thus, assuming a nonresonant £l r , the 
recreated torus to O(e) is located at 



J(9,t) = J r -d e Vb(9). 



(24) 
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B. Simplifying the control term 

The control fields f(6,t) and f a (0,t) obtained above 
can be used for studying the classical dissociation dynam- 
ics. However, a direct use of the control terms in quantum 
studies is subtle since the notion of action-angle variables 
does not exist except in the semiclassical limit. Thus, in 
order to implement the classical control terms for study- 
ing their effect on the quantum dissociation dynamics it 
is necessary to simplify the form of the control field. For- 
tunately, Huang et al. have already suggested .53] such a 
simplification and we briefly outline their approach. 

The control term, being periodic in 6 and t has a rich 
Fourier spectrum. However, only few of the Fourier com- 
ponents are dominant and the parameter 



Gk ± ,k 2 



\kiQ r + k 2 u>F\ 



(25) 



with -Ffc^fca being the coefficients in the double Fourier 
expansion of the control term / or f a is used to iden- 
tify those dominant modes. Note that this implies large 
amplitude Fk lt k 2 an d kiQ r + k 2 top ps i.e., the corre- 
sponding wavevector is close to being in resonance with 
the frequency vector u> of the integrable motion. Once 
identified, only the dominant Fourier modes are retained 
in the control term. Further simplification, required for 
quantum studies, is obtained by mapping a typical dom- 
inant term as 

F klj k 2 cos(fci# + k 2 Lu F t) — > A feli fc 2 cos(k 2 u) F t). (26) 

The coefficients Xki,k 2 are determined [53j by comparing 
the dominant Fourier mode amplitudes in the original 
control Hamiltonian with the corresponding amplitudes 
in the simplified control Hamiltonian 



H(J, 0- 1) + fx(x)X ku k 2 cos{k 2 ujFt). (27) 



If more than one dominant Fourier modes are present 
then they will appear as additional terms in Equationl2"Tl 
For all the results presented in the next few sections 
we have used the control Hamiltonian of the form given 
above. 



V. INFLUENCE OF THE CONTROL FIELD ON 
DISSOCIATION DYNAMICS 

We now present our results for the effect of local phase 
space barriers on the dissociation dynamics of the Morse 
state v = 10 for the two representative field frequencies 
and compare to the uncontrolled results summarized in 
Fig. [T] and Fig. [2] As we are interested in understanding 
the effect of creating cantori barriers on both the classical 
and quantum dynamics we also show, follo wing ea rlier 
studies on quantum transport through cantori[46, 64, 65], 
the time-averaged probability 



1 f T 

P v y = hm — / drKxv 



(28) 



of being in a state Xv' having started in the initial state 
X v . In this work T — 500r is a sufficiently long time 
for computing P v y . The classical analog of Eq. [28] is 
constructed by coarse-graining the actions i.e., the tra- 
jectory is considered to be in the action region J if it is lo- 
cated within a bin of width 0.5 centered about J. Reason- 
able variations of the bin width lead to qualitatively sim- 
ilar results and convergence can be easily checked. Such 
a coarse-graining procedure is appropriate for studying 
the classical-quantum correspondence of P v y . 



A. Off-resonant laser field 

Figure [3] summarizes our results for the off-resonant 
case with two different cantori barrier being rebuilt in the 
phase space. These cantori, corresponding to ujp/£l r = 
1+7 _1 (shown in red), and l + 7~ 2 (shown in green), are 
located at actions J r ~ 13.8 and 12.0 respectively. The 
modified phase spaces shown in Figs.[3fc) and (d) clearly 
show the reconstruction of the respective barriers as ev- 
idenced by the reduction of stochasticity and increased 
stickiness around the regular regions. These phase spaces 
should be compared to the one shown in Fig. [2b) and 
as anticipated the control field strength in Eq. [57] A2 is 
indeed smaller than the driving field strength Ai. Specif- 
ically, A 2 » 0.017 for the 1 + j- 1 barrier and A 2 ps 0.011 
for the 1 + 7~ 2 barrier with the dominant Fourier mode 
being (ki,k 2 ) = (3, —2) in both cases. In other words the 
control field is of the form A2 cos(2wf£) and, since A2 > 0, 
comes with a phase difference of it relative to the driving 
field. Interestingly, Wu et al. in an earlier work[23] have 
suggested precisely the same control field characteristics 
for suppressing chaos in the driven Morse system. How- 
ever, they were not clear about the mechanism for the 
suppression and this work yields the necessary insight in 
terms of the creation of local cantori barriers. 

For further insights into the role and efficiency of the 
cantori barriers towards controlling the dissociation dy- 
namics, in Fig. OUa) and (b) the classical and quantum 
time averaged-probabilities P u y defined in Eq. [28] are 
shown. Also shown in these figures are the approximate 
locations of the cantori as thin vertical lines at the cor- 
responding action values J = J r . The classical P u y in 
Fig. [3ta) show that the probabilities fall rapidly in the 
vicinity of the rebuilt cantori, especially in case of the 
1 + 7~ 2 cantorus. Consequently, dramatic reduction of 
the classical dissociation probability in both the cases 
can be seen (circles) from Figs. [He) and (f). It is pos- 
sible to investigate more detailed aspects of the classi- 
cal phase space transport across the cantori, as done for 
other systems (fril [65|. but we do not pursue them in this 
work. Moreover, it is known that the driven Morse oscil- 
lator dynamics near the separatrix can be analyzed from 
the perspective of a whisker-map for which Maitra and 
Heller have already provided a detailed classical-quantum 
correspondence of the transport across cantori[65j. 

A key issue that we are interested in this paper is 
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FIG. 3: (Color online) Effect of creating local barriers 
luf /£l(J) = 1 + 7 _1 (red, phase space shown in (d)) and 
ljf/Q(J) — 1 + 7~ 2 (green, phase space shown in (c)) on 
the classical (circles) and quantum (squares) dissociation dy- 
namics for luf = 0.0129 (off-resonance). The perturbatively 
created tori according to Eq. [M] are shown in the respec- 
tive phase spaces. Panels (a) and (b) show the classical and 
quantum time-averaged (T = 500r) cross probabilities (cf. 
Eq. I28J1 respectively with the initial state being the Morse 
state v = 10. The thin vertical lines show the expected lo- 
cation of the KAM barriers. In (e) and (f) the dissociation 
probabilities are shown as a function of time and correspond 
to (d) and (c) respectively. In the entire figure, filled and open 
(red, green) symbols are for the uncontrolled and controlled 
cases respectively. See text for discussion. 



whether the quantum dissociation is sensitive to the clas- 
sical phase space barriers being rebuilt. Note that all the 
quantum calculations performed herein have h = 1 and 
hence we are in the 'quantum regime'. Therefore, a pri- 
ori one might anticipate that quantum effects can over- 
ride or ignore the changes in the classical phase space. 
However, in this off-resonant case, we see from Figs.[3][e) 
and (f) that the quantum results (squares) exhibit clear 
reduction in the dissociation probabilities (See inset to 
Fig. [D^a) for the entire range of field frequency) . Anal- 
ogous to the classical case, the I + 7~ 2 cantorus is a 
stronger barrier to dissociation as seen by comparing 
Fig. [3Je) with Fig. [3Jf). The quantum time-averaged 
probability P v y is shown in Fig.[3]Jb) and exhibits the ex- 
pected suppression of probabilities for states lying around 
and beyond the location of the classical cantori. Com- 
paring the quantum P v y with the classical results shown 
in Fig.[3Ja) we make a few important observations. First, 
the finite probabilities for low lying Morse states (y < 4) 
seen in the quantum P v y are strictly zero in the classi- 
cal case. This is due to dynamical tunneling through the 
classical KAM barriers as proposed nearly two decades 



ago by Davis and Wyatt[I3|. Second, the quantum re- 
sults exhibit oscillations beyond the cantori barrier in 
contrast to the smooth classical decay. We suggest that 
this is a manifestation of what Maitra and Heller called 
'retunncling' in their study [6a] of the whisker map. Al- 
though the reconstructed cantori are perceived as com- 
plete barriers by the quantum system, some of the quan- 
tum states are able to tunnel efficiently across the cantori 
since Ti is large i.e., the quantum mechanism (enhance- 
ment due to tunneling) dominates the classical (suppres- 
sion due to cantorus ) mechanism. This might explain as 
to why the suppression of quantum dissociation probabil- 
ity due to the 1 + 7 -1 barrier is not significantly different 
from that due to the I + 7 -2 barrier in contrast to the 
classical results. 

Despite the comments made above, it is clear from 
Fig. [3] that the classical-quantum correspondence holds 
for local phase space control in the off-resonance case. 
We now discuss the on-resonant case wherein such a cor- 
respondence is not expected to hold. 



B. On- resonance laser field 

As mentioned in the previous section, with the primary 
driving field frequency value of u>f — 0.0178 the quantum 
system is in the two-photon resonant regime involving the 
Morse states v = 10, 12, and 14. This is reflected in the 
quantum dissociation probabilities shown in Fig. [UJd) as 
well as in the classical phase space as a large 2:1 nonlin- 
ear resonance zone (cf. Fig. [2Kb)) . Importance of this 
resonance, embedded in the chaotic region between the 
unperturbed 1 + 7 -1 and 2 + 7 _1 cantori, was noted by 
Brown and Wyatt[lj| as well as Dardi and Gray[l7| in an 
earlier work. Indeed, our computations (not shown here) 
indicate that the Wigner function of state v = 12 is lo- 
calized on the resonance with the Wigner functions asso- 
ciated with v = 10 and 14 straddling the resonance zone. 
Such a situation is tailor made for the manifestation of 
resonance-assisted tunneling in the system [47l IBlil. l57l. [5c| . 
Combined with the observation that the initial state of 
interest is located right around the 1 + 7" 1 cantorus (see 
Fig. H{c)), one expects significant competition between 
the quantum and classical mechanisms for dissociation. 
Consequently, the two-photon case provides a difficult 
challenge for the local phase space control method. 

In Fig.[¥]we summarize the results of our attempts to 
control the dissociation by creating the local barriers with 
cop/rir w v3 (red) and Wj?/f2 r = 1+7" 1 (green). We did 
not attempt to create the l+7~ 2 barrier since it would be 
located much below the action (J = 10.5) of the initial 
state in the classical phase space. Results for the two 
cases will be discussed separately in order to illustrate 
the interplay between classical and quantum dissociation 
mechanisms. Moreover, in the case of the ujf/^t ~ VS 
barrier complications arise in determining the simplified 
control Hamiltonian which requires additional discussion. 

Since Fig. [2f c) shows that the initial state is located 
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FIG. 4: (Color online) Results for the on-resonant case 
ujf = 0.0178, with notations as in Fig. [3] showing the ef- 
fect of creating local barriers ujf /fi(J) w %/3 (red) and 
ojf/£1(J) = 1 + 7 _1 (green). The respective phase spaces 
near the 2:1 resonance are shown in (c) and (d). Note that 
in this figure the results corresponding to ujf/0(J) ~ v3 are 
obtained by retaining the (3, —2) Fourier mode alone. See the 
discussion following Eg. 1291 and Fig. [5] for details. 



in the region corresponding to the 1 + 7 _1 cantorus we 
attempt to rebuild the cantorus and Fig.[4jd) shows that 
the attempt is successful. The control field in Eq. [57] is 
of the form A2 cos(2wpi)) corresponding to the dominant 
Fourier mode ^3,-2, with A2 ~ 0.008 and hence, as in the 
previous off-resonance case comes with a relative phase of 
7r. Robust barriers have been created in the phase space 
with a local control field strength A2 which is nearly four 
times weaker than the driving field strength. The time- 
averaged probabilities in Fig. [Ha) exhibit rapid decay 
in the vicinity of the recreated barrier and the classical 
dissociation probability, shown in Fig. [4fd), is reduced 
by nearly a factor of two. Surprisingly enough, Fig. H^e) 
shows that the quantum dissociation is enhanced by a 
small amount consistent with the behavior of the time- 
averaged probabilities shown in Fig. . The quantum 
result, in contrast to the off-resonance case, indicates that 
both classical and quantum mechanisms are at work in 
this instance. 

On the other hand, attempts to create the uif/^t ~ 
V3 barrier poses a problem, associated with the simpli- 
fication of the control term Eq. [27] which was not en- 
countered in the previous examples. Interestingly, in this 
case two dominant Fourier modes i*3 _ 2 and ^4,-2 are 
found with the corresponding values G3-2 ~ 0.029 and 
G 4i _ 2 ~ 0.024 (cf. Eq. 123]), Taking into account only 
the marginally dominant (3, —2) mode Fig. [JJc) shows 
that the desired local barrier is created. Consequently, 



Fig. Q|f) shows that the classical dissociation probabil- 
ity is, as in the case of the 1 + 7 _1 barrier, reduced by 
a factor of two. Again one observes that the quantum 
counterpart behaves in an opposite manner i.e., the dis- 
sociation is slightly enhanced. However, given that the 
Fourier mode (4, —2) is nearly as dominant as the (3, —2) 
mode, it seems reasonable to use an effective A2 in the 
simplified control term of Eq. [371 as 



A, 



3,-2 



-F4.-5 



v 3 (j r ) v 4 (j r y 



(29) 



Such a procedure yields A2 ~ —0.015 and thus the con- 
trol field, still less intense than the primary field, comes 
with a relative phase of zero. Nevertheless, such an at- 
tempt fails as seen by inspecting the phase space shown 
in Fig. Etc) which exhibits increased stochasticity and, 
expectedly, leads to enhanced classical dissociation ob- 
served in Fig. E{d). Now, however, Fig. (5je) shows that 
the quantum dissociation is suppressed appreciably and 
hence it would seem as if the quantum dynamics feels 
the barrier where there is none! Arguments invoking the 
large, finite value of fi and dynamical localization can be 
safely ignored since the phase space in Fig. (5jc) does not 
show any appreciable stickiness near the apparent bar- 
rier. The resolution to such an unexpected result comes 
from a closer inspection of the phase space in Fig. EJc). 
One can clearly see that the attempt to create the lo- 
cal barrier has resulted in a severe perturbation of the 
2:1 nonlinear resonance. As the additional Fourier mode 
(4, —2) implies 2il r u>p, the observed perturbation can 
be tracked to the specific mode as long as it comes with 
an opposite phase. This will be established in the follow- 
ing subsection (cf. Eq. 131)]) . Since this specific resonance 
is key to the two-photon process, it should not be sur- 
prising that the quantum dissociation is suppressed. 

It is crucial to note that failure to create the barrier 
of interest occurs only when the control term is mapped 
to the simplified form as in Eq. [27] using the effective 
value for A2 shown above. This is confirmed by inspecting 
the phase spaces shown in Fig. [Sfa) , (b) associated with 
the full 0(e 2 ) control Hamiltonian (see Eq. [22]) and the 
approximate Hamiltonian 

H c {j,e,t) » H(j,e,t) 

+ F « -2cos(n0 - 2u F t), (30) 

n=3A 

obtained by retaining only the dominant Fourier modes 
in the control term. The recreated barriers in the phase 
space can now be clearly seen and Fig. [H^d) shows that 
the classical dissociation computed using the Hamilto- 
nian in Eq. [3D] is indeed suppressed. Clearly, the oppos- 
ing classical and quantum results in this subsection, with 
the associated phase space structures, point to the im- 
portance of the quantum dissociation mechanism in the 
on-resonant case. In what follows we show that these ob- 
servations can be rationalized based on the phenomenon 
of resonance-assisted tunneling. 
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FIG. 5: (Color online) Phase spaces corresponding to creat- 
ing wp/f2(J) « \/3 barrier in the on-resonant case with (a) 
Full leading order control term (cf. Eq. I22[) (b) Retaining two 
dominant Fourier modes (cf. Eq. I30[) and (c) Simplified con- 
trol term (cf. Eq. I27[) using both Fourier modes and effective 
strength A2 obtained from Eq. I29I In (a) the size of the black 
square represents h — 1. The rebuilt torus can be seen in (a) 
and (b), highlighted in green, but not in (c). Note the per- 
turbation of the 2:1 resonance in case (c) indicated by a red 
arrow. In (d) and (e) the classical (corresponding to the figure 
(b) only) and quantum dissociation probabilities respectively 
are shown for the various cases. Curves marked by red arrow 
correspond to the phase space (c). 



1. Dissociation via resonance- assisted tunneling 

In order to confirm the above arguments and to gain 
a better understanding of the results shown in Fig. [4] 
and Fig. O for the dissociation of the state v = 10 we 
focus on the role of 2:1 resonance within the paradigm 
of resonance-assisted tunneling. Following the earlier 
works [5?1 [57l . [58l | . the motion in the vicinity of a r:s 
resonance is analyzed by applying secular perturbation 
theory on the uncontrolled Hamiltonian in Eq . [71 and for 
details we refer the reader to the original work [561]. First, 
a canonical transformation to the appropriate slow angle 
9 — > cf) = 9 — Q r:s t, J — ► J is made resulting in the new 
Hamiltonian 

H(J,cf>,t) = H (J) - Sl r:s J + V(J,<i>,t), (31) 

where V(J, (j), t) = V(J, </> + Q r:s t, t) and we have denoted 
V(J,9,t) = -ev{J,6;t) (Cf. Eq. EJ. We now expand 
Hq in Eq. [3T] about the resonant action J r:s to obtain 
the zeroth-order Hamiltonian in the vicinity of r:s of the 
form 



Hq[J) — Ho[J r: , 



2m r - 



(32) 



with AJ = J—J r:s and m r:s = —2Dq/lj'q. Since, 4> varies 
slowly near r:s, V(J, (j), t) is replaced by its average over 
r field periods 

Kv(J,0) = — I V(J,<t>,i)dt 



Ai^K(J)cosH). (33) 



Ignoring the higher harmonics in the above expression 
and neglecting the action dependence of the Fourier co- 
efficients V n we obtain an effective (integrable) pendulum 
Hamiltonian 



Heff(J,<f>) 



1 



2rh r 



-(AJ) 2 + 2V; :s (J r:s )cos(r0), (34) 



to describe the dynamics near the r:s resonance with 
rhr-.s — —Tn r . s and V r:s (J) — XidiV\(J)/2. 

Specializing the above result to the observed 2:1 res- 
onance, the resonant action is determined to be J2:i — 
12.6, thus confirming the participation of the Morse state 
v = 12. Using the zeroth-order Hamiltonian in Eq. I3~4l 
one finds, with J = 10.5 (quantum state v = 10) and 
J' = 14.5 (quantum state v' = 14), that 



\Er - E, 



1 



2m 2: i 
3.2 x 10" 4 



(j- j'-2j 2:1 ; 



(35) 



i.e., the states v — 10 and v = 14 are nearly sym- 
metrical with respect to the state v — 12 localized on 
the 2:1 resonance (47l [56| . Therefore, the nonzero cou- 
pling V<x;\ will efficiently connect the states v = 10 
and v — 14. Moreover, it is possible to estimate the 
strength of the resonance for the given parameters as 
V2:i(-/2:i) ~ 0.01464. It is crucial to note that the the ef- 
fective control field coupling strength from Eq.[^H]in case 
of ujp/fl r fa V3 i.e., relevant to the phase space shown 
in Fig. |Sfc) satisfies 



A2 ~ — V2:i(J2:l)- 



(36) 



Thus the control field with a zero relative phase tends 
to negate the effect of the 2:1 resonance generated by 
the primary driving field and this can be clearly seen 
in Fig. EJc). We also remark here that the substantial 
quantum dissociation probability seen in Fig. [5Je), de- 
spite the strong perturbation of the 2:1 resonance, is due 
to the higher harmonics which have been neglected in 
the present analysis. This confirms our suspicion that 
the decay of state v = 10 in the two-photon regime is 
dominated by resonance-assisted tunneling. 



C. Are the rebuilt barriers local? 

Up until now most of our results and discussions have 
focused on a specific initial Morse state. The barriers 
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FIG. 6: (Color online) Effect of the local control term on 
the dissociation probability of nearby states, specified by v, 
at t — 500r. (a) and (b), with off-resonance uif = 0.0129, 
correspond to the creation of the 1 + 7 -1 and 1 + y~ 2 bar- 
riers respectively. Panels (c) and (d), with on-resonance 
lof = 0.0178, correspond to the creation of the 1 + j^ 1 and 
\/3 barriers respectively, with fixed field strength as in the 
previous figures. The uncontrolled results are shown as filled 
symbols, classical (circles) and quantum (squares), while the 
controlled results are shown as corresponding open symbols. 
The lines are drawn as a guide to the eye. 



were created locally in phase space to influence the dis- 
sociation dynamics of the state v = 10. From a control 
point of view, it is of some interest to examine the effect 
of such barriers on the dynamics of other states, espe- 
cially states that are in the vicinity of state v = 10. 
In other words, how local are these barriers? Towards 
this end, in Fig. [5] we show the influence of the local 
barriers on the dissociation probabilities of other nearby 
Morse states with v > 10. At this stage it is usefull to 
recall the results shown in Fig. [1] with the essential dif- 
ferences between the off-resonant and on-resonant cases. 
The effect of the 1 + 7 -1 and 1 + 7~ 2 barriers in the 
off-resonant case, located around J r ss 13.8 and 12.0, are 
shown in Figs. [Sfa) and (b) respectively. These figures 
confirm that to a large extent the barriers are indeed local 
i.e., dissociation is suppressed to varying extent for ini- 
tial states lying below the barrier. For states lying above 
the barrier, dissociation is either enhanced (mostly in the 
quantum case) or slightly reduced. 

On the other hand the quantum results for the on- 
resonant case, shown in Figs.^c) and (d) for the I + 7 -1 
(J r ss 9.9) and \/3 (J r ~ 10.9) barriers (correspond- 
ing to Fig. \M,,c) i.e., with only the (3, —2) Fourier mode 
included) respectively, are far less straightforward to in- 
terpret. This is not entirely surprising since, as shown 



in the last section, resonance-assisted tunneling plays an 
important role and overrides the suppression due to the 
local barriers. In particular, the dissociation probabili- 
ties for states v = 10 and v = 14 increase very slightly. 
However, for the states v = 11 and v = 13, not in- 
volved in the resonance-assisted tunneling process, one 
observes substantially reduced dissociation despite being 
located far above the local barriers. We suspect that 
this is due to the increased stickiness around the 2:1 res- 
onance region observed in the controlled phase spaces 
shown in Figs. BJc) an d (°0- Note that this is corrobo- 
rated by the observation that the concerned states also 
exhibit reduced classical dissociation as seen in Figs. [He) 
and (d). Further confirmation comes from our calcula- 
tions which show an opposite quantum trend to that of 
Fig. Od) upon inclusion of the (4, —2) Fourier mode as 
well resulting in the phase space shown in Fig. [H^c). It is 
possible to implicate, albeit indirectly, the local barriers 
with the observed suppression since creating the barriers 
leads to a more regular phase space, slightly increased 
size of the 2:1 resonance region, and therefore increased 
stickiness. Nevertheless, comparing the off-resonant and 
on-resonant cases shown in Fig. [51 it is evident that the 
effects of creating local phase space barriers can be far 
more subtle to interpret in the latter case. 



VI. CONCLUSIONS 

To summarize, this work demonstrates that it is pos- 
sible to control the dissociation dynamics of a driven 
Morse oscillator by creating local phase space barriers. 
A clear understanding of the effect of cantori on both 
classical and quantum dissociation dynamics is obtained 
(see inset in Fig.^a) for example). This work also high- 
lights the essential difference between off-resonant and 
on-resonant dynamics with resonance-assisted tunneling 
playing a prime role in the latter case. Although local 
phase space barriers are very efficient in reducing the dis- 
sociation in the off-resonant case, the results in Fig.[4]and 
Fig. [5] suggest that controlling the quantum dissociation 
in the on-resonant regime can be achieved by using con- 
trol fields which selectively perturb the appropriate non- 
linear field-matter resonance. Similar observations have 
been made earlier in a different context [58| and further 
studies are required to establish a general criteria for con- 
trolling the multiphoton processes from the viewpoint of 
local modification of the phase space structures. 

Several questions, however, remain to be addressed and 
we mention a few important ones. First, there is the issue 
of the effectiveness of the local barriers in systems with 
more than two degrees of freedom since the invariant tori 
do not have the correct dimensionality to partition the 
phase space. However, there are reasons to hope that 
even for systems with higher degrees of freedom the re- 
built tori can act as barriers for short times [6^]. A care- 
ful study of the classical and quantum transport with 
and without the control fields is required in this instance. 
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Apart from polyatomic molecules, this is also important 
while considering the infrared multiphoton dissociation 
dynamics of a Morse oscillator by explicitly taking into 
account the rotations [o7l [68j . Second, one would like 
to extend the approach to systems under the influence 
of more general time-dependent fields, for e.g., chirped 
fields. For slow chirping this should be feasible as one can 
then utilize the concept of adiabatic Floquet theory (69|. 
In case of arbitrary time-dependent fields the correct ap- 
proach is not clear at the present moment. Third, the 
method used in the present and earlier works is depen- 
dent on our ability to express the Hamiltonian in terms 
of appropriate action-angle variables. This might pose 
some limitations which, as seen in the present work, can 
be more severe in terms of implementing the control into 



the quantum dynamics. Finally, for systems with small 
effective h, the resonance-assisted tunneling mechanism 
will be replaced by chaos- assisted tunneling [7(| and it 
would be interesting to study the influence of the cantori 
barriers in this context. These issues are the focus of 
ongoing work in our group. 
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